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We present deep images of dust continuum emission at 450, 800, and 850 //m 
of the dark cloud LDN 1689N which harbors the low-mass young stellar objects 
(YSOs) IRAS 16293-2422A and B (I16293A and I16293B) and the cold prestel- 
lar object I16293E. Toward the positions of I16293A and E we also obtained 
spectra of CO-isotopomers and deep submillimeter observations of chemically re- 
lated molecules with high critical densities (HCO+, Hi3C0+, DCO+, H2O, HDO, 
H2D+). To I16293A we report the detection of the HDO Iqi-Ooo and H2O Iiq- 
loi ground-state transitions as broad self-reversed emission profiles with narrow 
absorption, and a tentative detection of H2D+ lio-lii- To I16293E we detect 
weak emission of subthermally excited HDO Iqi-Ooo- Based on this set of sub- 
millimeter continuum and line data we model the envelopes around II 6293 A and 
E. The density and velocity structure of I16293A is fit by an inside-out collapse 
model, yielding a sound speed of a=0.7 km s^^, an age of t=(0.6-2.5) x 10"^ yr, 
and a mass of 6.1 M0. The density in the envelope of I16293E is fit by a radial 
power law with index —1.0 ± 0.2, a mass of 4.4 Mq, and a constant temperature 
of 16 K. These respective models are used to study the chemistry of the envelopes 
of these pre- and protostellar objects. 

We made a large, fully sampled CO J=2-l map of LDN 1689N which clearly 
shows the two outflows from I16293A and B, and the interaction of one of the flows 
with I16293E. An outflow from I16293E as reported elsewhere is not confirmed. 
Instead, we find that the motions around I16293E identified from small maps are 
part of a larger scale fossil fiow from I16293B. Modeling of the I16293A outfiow 
shows that the broad HDO, water ground-state, and CO J=6-5 and 7-6 emission 
lines originate in this flow, while the HDO and H2O line cores originate in the 
envelope. The narrow absorption feature in the ground-state water line is due to 
cold gas in the outer envelope. The derived H2O abundance is 3 x 10~^ in the cold 
regions of the envelope of I16293A (Tkin < 14 K), 2 x 10~^ is warmer regions of 
the envelope (> 14 K), and 10~* in the outflow. The HDO abundance is constant 
at a few xlO"^*^ throughout the envelopes of I16293A and E. Because the derived 
H2O and HDO abundances in the two objects can be understood through shock 
chemistry in the outflow and ion-molecule chemistry in the envelopes, we argue 
that both objects are related in chemical evolution. The [HD0]/[H20] abundance 
ratio in the warm inner envelope of I16293A of a few times 10~^ is comparable 
to that measured in comets. This supports the idea that the [HD0]/[H20] ratio 
is determined in the cold prestellar core phase and conserved throughout the 
formation process of low-mass stars and planets. 

Subject headings: stars: formation — ISM: clouds — ISM: jets and outflows — 
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ISM: molecules — ISM: individual (IRAS16293-2422) — astrochemistry 



1. Introduction 

IRAS 16293—2422 is one of the best studied low-mass young stellar objects (YSOs) 
(Blake et al. 1994; van Dishoeck et al. 1995; Ceccarelli et al. 2000). It is deeply embedded 
in the LDN 1689N cloud in Ophiuchus (distance 160 pc) and classified as an extreme Class 
I/Class YSO (Lada 1991). It is a low-luminosity (27 Lq, Walker et al. 1986) binary system 
(hereafter, 116293A and 116293B) with a separation of 840 AU (5") and individual stellar 
masses of 0.5 M© (Mundy et al. 1992). Two powerful outflows emanate from the system. At 
least three physically and chemically different regions can be recognised in a 20" (3000 AU) 
beam (van Dishoeck et al. 1995): (i) a cold, relatively low-density outer envelope (Tkin ~10- 
20 K, n(H2) ~ 10^-10^ cm~^) which marks the transition into the extended parental cloud 
LDN 1689N; (ii) a warmer circumbinary envelope of size ~ 2000 AU (Tkin — 40 K, ti(H2) ~ 
10^-10^ cm-3); and (iii) a warm, dense core of about 500-1500 AU (Tkin > 80 n(H2) ~ 10^ 
cm~^) which traces the interaction of the outflow and the stellar radiation with the inner 
part of the circumbinary envelope. This picture emerged from excitation analysis of species 
tracing these distinct regions and conditions, assuming uniform temperatures and densities 
for each. Continuous descriptions of the density and temperature as function of radius in the 
collapsing envelope were developed by CeccareUi et al. (2000) and Schoier et al. (2002) on 
scales from 30 to 5000 AU. In this Paper, we present new observations and modeling, that 
allows us to independently derive the density and temperature structure in the extended 
envelope around the (pre-) protostellar cores through continuum photometry and imaging 
of the dust, and heterodyne spectroscopy of the molecular gas. The resulting density and 
temperature distribution ranges from the warm inner envelope region (100 AU) to the coldest 
(12 K) outer envelope regions (7300 AU) and serves as a basis for excitation and abundance 
analysis of the line data. 

Another component of the IRAS 16293—2422 system is located 80" east and 60" south of 
116293A. This smafl clump, 116293E, has strong DCO+, NH3, and NH2D emission (Wootten 
& Loren 1987; Mundy et al. 1990; Shah & Wootten 2001) connected to the 116293A core, 
and signiflcant submillimeter continuum emission from dust (Sandell 1994). Single-dish 
observations of NH3 in a 40" beam indicate a temperature Tkin = 12 K (Mizuno et al. 1990). 
Millimeter-interferometric data, sampling smaller scales and less sensitive to the warmer 
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extended envelope, yield Tkin ~ 8 K and n(H2) ~ 1 x 10^ cm~^ if homogeneous conditions 
are assumed (Shah & Woottcn 2001). Because no far-infrared (FIR) or submillimeter point 
source has been found embedded in I16293E, it is classified as a pre-(proto)stellar core. In 
this Paper, we present submillimeter images that constrain the density structure of I16293E. 

It is difficult to study chemical evolution in the early protostellar stages, because during 
coUpase the temperature stays low at ~ 10 K while the density increases to n(H2) > 10^ 
cm~^. Under these conditions the line emission of many molecules is dominated by a warmer 
outer envelope around the prestellar core. Also, many molecules in the cold core will condense 
on dust grains. Eventually, only H2, H^, and their isotopomers HD and H2D+ will remain 
in the gas phase. This leads to a significant enhancement of H2D"'" through the deuterium 
exchange reaction 

H+ + HD ^ H2D+ + H2 + AE, (1) 

since the backward reaction becomes neghgible at low temperatures (Smith et al. 1982; 
Herbst 1982). H2D+ is thought to play a pivotal role in the deuteration of molecules (Watson 
et al. 1976; Millar et al. 1989) and the detection of H2D+ in the Class YSO NGC 1333 
IRAS 4A (Stark et al. 1999) is an important confirmation of the cold gas-phase chemical 

networks. The H2D"'" enhancement is reflected in the high abundance ratios that are generally 
observed of, e.g., [DCO+]/[HCO+], [NH2D]/[NH3], [DCN]/[HCN], and [N2D+]/[N2H+] in 
cold prestellar cores and Class YSOs (e.g., Cuelin et al. 1982; Olberg et al. 1985; Wootten 
1987; Butner et al. 1995; Williams et al. 1998; Bergin et al. 2002). Models predict that in the 
collapse stage the abundances of many deuterated radicals and molecules increase sharply 
after their non-deuterated versions start to get heavily depleted on dust grains, then reach a 
peak and start to decrease as the deuterated molecules too condense on the dust (e.g.. Brown 
& Millar 1989; Rodgers & Millar 1996). Reactions with atomic deuterium on the grain surface 
may further enhance the abundance of deuterated species in the ice (Tielens 1989; Brown 
& Millar 1989). After the formation of a YSO, the condensed molecules may evaporate 
from grains, resulting in a temporary increase in the gas-phase deuteration fraction. The 
deuteration may decline after ~ 10"^ years (Rodgers & Millar 1996) when higher temperature 
chemistry becomes effective. In addition, in such hot (Tkin — 100-200 K) regions the reversal 
of reaction (1) becomes dominant and very little new fractionation is expected to occur. We 
report in this Paper observations of several key deuterated species in I16293A and I16293E. 

Water is an important species in the early stages of star formation. H2O is the most 
important hydride molecule in oxygen chemistry and models predict that H2O plays a dom- 
inant role in the energy balance and chemical evolution during star formation. Water will 
have the largest abundance and excitation contrast between the protostellar source and the 
surrounding cloud (e.g., Ceccarelli et al. 1996; Doty and Neufeld 1997). Measurements of the 
rotational and ro- vibrational transitions of H2O in star- forming regions have recently become 
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available from the Infared Space Observatory (ISO; e.g., Helmich et al. 1996a; Ceccarelli et 
al. 2000) and the Submihimeter Wave Astronomical Satelhte (SWAS; e.g., Ashby et al. 2000; 
Snell et al. 2000; Neufcld ctal. 2000). In particular the heterodyne velocity-resolved measure- 
ments of the ground-state transition of ortho-H20 with SWAS allow a direct determination 
of the water abundance throughout the envelope. The HDO Ioi^Oqo ground-state transition 
at 465 GHz {E^ = 23.2 K, rzcrit > 10^ cm'^ for Tkin < 50 K, Green 1989) is an excellent 
tracer of high density gas at low temperatures where the abundance of HDO relative to water 
is enhanced. Ground-based observations of the HDO ground-state transition are still very 
sparse. To date, this transition has only been observed in the high-mass star-forming regions 
Orion-KL (Schulz et al. 1991), W3(0H), and W3(H20) (Helmich et al. 1996b). The latter 
study indicates that the [HD0]/[H20] abundance ratio is comparable to that found in hot 
cores and is not a sensitive indicator of the evolutionary stage in high-mass star formation. 
This could mean either that the W3 cloud always stayed warm or that the [HD0]/[H20] ratio 
rctrusn to thermal equilibrium faster than ratios like [DCN]/[HCN] (Helmich et al. 1996b). 
For a low-mass Class YSO like 116293A the situation may be different, because the time 
scale to reach steady-state may be much longer. In this Paper, we present observations of 
HDO and H2O of I16293A and E. 

In this paper we present a detailed study of the physical and chemical structure of 
the low-mass star-forming cloud LDN 1689N with particular emphasis on the H2O/HDO 
chemistry and deuterium fractionation. We present extensive submillimeter continuum pho- 
tometry and maps reveahng the warm and cold dust around I16293A, B and I16293E (§2.1). 
We report the detection of the ground-state lines of (deuterated) water and a tentative detec- 
tion of H2D+ toward I16293A as well as weak HDO emission in the prcstcllar core I16293E 
(§2.2). These data arc combined with pointed observations of CO and HCO+ isotopomcrs to 
determine the temperature and density structure throughout the envelope (§3). We address 
the two outflows and their origin in §4. Our deep observations of the HDO and H2O ground- 
state transitions are used to study the [HD0]/[H20] ratio throughout the envelope (§5.1). 
A simple chemistry network is used to model the deuterium chemistry and in particular the 
water fractionation (§5.2). Finally, we discuss the evolutionary difference between I16293A, 
I16293B, and I16293E (§6), and summarize the conclusions in §7. 

2. Observations and Data Reduction 

2.1. Submillimeter Continuum 

All continuum observations were made with the James Clerk Maxwell Telescope (JCMT), 
using the single channel bolometer UKT14 and the versatile SubmiUimetre Common-User 
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Bolometer Array (SCUBA) under good weather conditions. The UKT14 measurements were 
done during the winters of 1991-1993. This instrument was for a long time the common- 
user bolometer system and is described in Duncan et al. (1990). It was replaced in 1997 by 
SCUBA which was designed both for photometry and fast mapping (Holland et al. 1998). 

The photometry of I16293E was done with UKT14 during Spring 1993. The chop throw 

was 90 " in Declination and the calibrators were Uranus and the nearby secondary calibration 
standard I16293A (Sandcll 1994). The results are given in Table 1. UKT14 maps were 
obtained at 1.1 mm, 800 /xm, 750 /im, 450 /xm, and 350 /im during several observing runs in 
1991 and 1992. Most of these maps have poor scnsivitity and were mainly used to derive the 
submillimctcr position of I16293A. However, even with the somewhat limited signal-to-noise 
(S/N) these provide a better estimate of the background subtracted flux of I16293A than 
single position photometry. 

At 800 iim we attempted to go deep and completely map the dust cloud surrounding 
the bright submillimeter source. Because the dust emission is very extended (see the 850 
lim SCUBA map in Fig. 1) we combined six maps of the region, centered on I16293A or 
I16293E. Each map was bound by pointing observations of the nearby blazars 1514—241 and 
1730—130. Typical pointing drifts were less than 2". These errors were removed in the data 
reduction process assuming that the pointing drifts linearly in Azimuth and Elevation as a 
function of time. All maps were made on-the-fly, scanned in Azimuth with a cell size of 4" 
and a chop throw of 40 " in the scan direction. With a single pixel instrument it is difficult 
to map such an extended region. The final map was therefore reduced with the Dual Beam 
Maximum Entropy (DBMEM) algorithm written by John Richer (Richer 1992) for JCMT. 
The DBMEM map was basehne subtracted and recalibrated by using the integrated intensity 
in a 60" x 60" area around I16293A from maps reduced with N0D2 (Haslam 1974; Emerson 
et al. 1979). The 800/im-map is shown in Fig. 1 (top panel). From comparison with 
SCUBA maps (see below), it is clear that the 800/xm-map has failed to recover the faint, 
relatively uniform, low-level emission surrounding both cloud cores. This does not affect 
the morphology of the map, but does have a large effect on integrated intensities of the 
surrounding LDN 1689N cloud. The low-level emission is missing because the 800 /xm-maps 
were taken with the same rather short chop throw, and the individual maps have a noise 
level comparable or higher to that of the low-level extended emission. It is also well known 
that the NOD2-algorithm is not very sensitive to relatively uniform extended emission on 
spatial scales of several times the chop throw (Emerson et al. 1979). 

Additional maps were obtained with SCUBA during Spring 1997 as part of commis- 
sioning jiggle- and scan-map observing modes. We have complemented these maps with 
additional maps obtained from the JCMT archive in the time period 1997-1998. All maps 
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were taken under excellent sky conditions (atmospheric opacities T225 ghz < 0.04), and all 
the maps in the archive were taken with a 120 " chop throw, typically in Azimuth. Most of 
the maps of I16293E were also obtained with a 120" chop. 

The pointing was checked before and after each map using the same blazars as for 
our UKT14 observations. Some of the archive maps, however, were taken without pointing 

observations. These were aligned with the pointing corrected average before adding them 
into the final data set. Calibration and beam characterization is based on beam maps of 
Uranus obtained with the same chop throw. The Uranus maps give a Half Power Beam 
Width of 14.5" and 7.8" for 850 /xm and 450 /xm, respectively. 

Both the 116293A, B and I16293E cloud cores are rather extended and no map is com- 
pletely free of emission in the off positions. Any bolometer that showed evidence for or was 
suspected of being contaminated by the chop was therefore blanked in the data reduction. 
Although this was done as carefully as possible, it is clear that the outskirts of our final 
mosaiced image (Fig. 1) may still be affected, especially since the dust emission extends 
beyond the area we mapped to the north and northeast of both I16293A and I16293E. The 
cloud core surrounding I16293E also extends further east than is covered by our maps. 

The basic data reduction was done using the SCUBA software reduction package SURF 
(Jenness 1998) as explained in the SCUBA Map Reduction Cookbook (Sandell 2001). In 
total we used 16 data sets for the 850 /im map and 13 for the 450 /im map. In the final 
coadd we adjust the pointing in each map (shift and add) to ensure that the final map is 
not broadened by small pointing errors. The rms noise levels in the final maps are difficult 
to estimate, because there are no emission-free regions in the maps, but we estimate them 
to be below 20 mJy beam~^ at 850 /im and below 80 mJy beam~^ at 450 /im. 

The final maps were converted to FITS and exported to MIRIAD (Sault 1995) for 
further analysis. In order to correct for the error beam contribution, especially at 450 /xm, 
we deconvolved the maps using clean and spherically symmetric model beams derived from 
beam maps of Uranus observed in stable night time conditions. At 850 /xm we use HPBWs of 
14.5", 55", and 150", with amplitudes of 0.985, 0.014 and 0.001, respectively. At 450 fim the 
HPBWs are 7.8", 34", and 140", with amplitudes of 0.964, 0.029, and 0.007, respectively 
The final maps were restored with a 14 " beam at 850 and an 8 " beam at 450 //m. The peak 
fiux densities in the restored maps are 16.0 Jy beam~^ and 76.1 Jy beam~^ toward I16293A 
at 850 and 450 //m, respectively, while the peak fiuxes of I16293E are 1.4 Jy beam~^ and 
4.35 Jy beam~^ at 850 and 450 //m, respectively. 
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2.2. Moleculcir line observations 

Spectra of CO J=2-l (230.538000 GHz) and isotopomcrs, as well as the J=3-2 transi- 
tions of DCO+ (216.112605 GHz), HCO+ (267.557619 GHz), and Hi3C0+ (260.255478 GHz) 
toward I16293A and I16293E were retrieved and reduced from a search of all released data 
in the JCMT archive at the Canadian Astronomy Data Center taken with the SIS receiver 
RxA2 and its successor RxAS. Deep observations of the HDO 3i2-2i2 (225.896720 GHz) and 
2ii-2i2 (241.561550 GHz) lines were made with the receiver RxA3 in 2001 April. We also 
acquired a large, fully sampled map in the CO 2-1 line which is centered on I16293E and 
covers the outflows from I16293A and B. The map was made on-the-fly with a bandwidth of 
125 MHz, a sampling of 5" in the scan direction, and a cross-scan step size of 10". The final 
map was made from a series of sub-maps scanned either in RA or Dec. The integration time 
per point was 5 seconds, but in poor weather conditions we often made multiple coverages. 
All observations were done in position switch mode with the reference position 600 " or 800 " 
west of the submiUimeter position of I16293E. At the same time we have also obtained spec- 
tra toward I16293A, which is a pointing and spectral line standard. The relative calibration 
of the spectra of the CO 2-1 map is accurate within 5%, because we overlapped the sub- 
maps. In 1999 April we obtained additional deep spectra of CO and -"^^CO 2-1 (220.398677 
GHz) at the peak positions in the outflow lobes in order to get more accurate estimates of 
the ^^CO optical depth in the high velocity gas. Most of the spectra were obtained with a 
bandwidth of 125 MHz, corresponding to a velocity resolution of 0.1km s~^. The spectra 
are calibrated in Tub (I16293A) or (I16293E) and shown in Figs. 2 and 3. The results 
of gaussian fits to the line profiles are summarised in Table 2. For some of the hues with 
symmetric profiles (e.g., HDO) the emission and the absorption components were fitted sep- 
arately. For most of the complex line profiles we only list the integrated emission and the 
velocity of the maximum self- absorption. Note that the phase-lock instability of RxA2 and 
RxA3 may cause an instrumental broadening up to 0.5 km s~^. 

The N2H+ J=4-3 (372.672509 GHz), H2D+ lio-ln (372.421340 GHz), DCO+ 5-4 
(360.169881 GHz), HCO+ 4-3 (356.734288 GHz), and }1^^C0+ 4-3 (346.998540 GHz) ro- 
tational transitions were observed in 2000 April with the JCMT using the dual channel 
SIS receiver RxB3. The receiver was operated in single side band (SSB) mode where the 
image side band is rejected and terminates in a cold load. The main beam efficiency at 
these frequencies is about 60%, while the beam width is 14-13 "for 345-372 GHz. The lines 
were observed in beamswitching mode switching over 180". The N2H+ and H2D+ lines were 
observed simultaneously with relatively good atmospheric transparency (r225 ghz — 0.05), 
but with only 30 minutes of useful integration time above an elevation of 40 degrees. RxB3 
was changed to tunerless mixers in 1999 and the performance near the upper band edge has 
decreased dramatically with respect to Stark et al. (1999). Nevertheless we have tentatively 
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detected the H2D+ ground-state line (Fig. 2). 

The HDO J=loi-Ooo ground-state transition (464.924520 GHz) was observed in 1998 July 
at the JCMT during a night of excellent submiUimeter transparency with a zenith optical 
depth at 225 GHz below 0.05. The dual channel, dual band 460/660 GHz SIS receiver RxW 
was used in SSB mode with the Digital Autocorrelator Spectrometer (DAS) backend. The 
DAS was split in two parts of 125 MHz with a spectral resolution of 189 kHz (=0.12 km s^^ 
at 465 GHz). The HDO line was observed in the upper side band, and SSB system temper- 
atures including atmospheric losses were Tgys ^ 1000 K. The beam size is about 11" at 465 
GHz and the main beam efficiency was measured to be about 50%. We observed the HDO 
line toward I16293A in position switch mode with the same reference as above. We also 
observed the hue in beam switch mode, switching over ±180" in Azimuth. These spectra 
show the same profile as the position switched ones, but only the beam switched spectra 
clearly reveal the level of the continuum emission from the dust heated by the YSO. All 
spectra were coadded and corrected for the continuum level. In addition, deep spectra of the 
HDO line toward I16293E were obtained in position switch mode. The data were calibrated 
similarly to the RxA2 and A3 observations. 

The H2O lio-loi ground-state transition (556.936002 GHz) was observed in 1999 August 
with SWAS (Melnick et al. 2000). The receiver consisted of a Schottky diode mixer resulting 
in a double sideband system temperature of about 2200 K. The back-end was an Acousto 
Optical Spectrometer with a 1.4 GHz bandwidth and a velocity resolution of < 1 km s~^. 
The beam width is about 4'. The H2O line was observed in position switch mode, switching 
to an emission free reference position 1.5° away. The total (on+off) integration time was 
23.5 hrs, yielding a noiose level of Tl(rms) = 0.01 K. 

The CO J=6-5 (691.473076 GHz) and 7-6 (806.651806 GHz) transitions toward I16293A 
have been taken with the JCMT using the RxW receiver in 1999 August and with the 
MPIfR/SRON 800 GHz SIS receiver in 2000 April, respectively, under good atmospheric 
conditions. The beam parameters at these frequencies were determined from observations of 
Mars. The main beam efficiencies are about 0.3 and 0.24, and the beam sizes (FWHM) are 
8" and 6" at 691 GHz and 806 GHz, respectively. 

3. Results and Analysis 

3.1. SubmiUimeter Continuum: Morphology 

IRAS 16293—2422 was found to be a binary system in 6 cm and 2 cm VLA observations 
by Wootten (1989). The two components are separated by 5" (840 AU) with a position 
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angle of PA — 135°. High resolution millimeter observations (Mundy et al. 1992; Looney 
et al. 2000) also resolve the IRAS source into a binary system, in good agreement with the 
VLA data. The southern component (A or MMl) is more extended in millimeter continuum 
and is associated with dense gas, high velocity emission and H2O masers (Mundy ct al. 1992; 
Choi et al. 1999; Looney ct al. 2000) suggesting that it is the more active component of the 
binary system. The northern component, B or MM2, is very compact both at centimeter 
and millimeter wavelengths, and has a surprisingly steep spectrum, a — 2.3 ± 0.3 (Mundy 
et al. 1992). It is brighter than A at both 1.3 cm and 2.7 mm (Estelella et al. 1991; Mundy 
et al. 1992; Schoier et al. 2003). 

We resolve IRAS 16293-2422 in all maps with a FWHM ~ 9.6" x 5.7" at 450 /xm with 
PA = 152°, and find it slightly more extended at 800 and 850 pim with PA = 145° and 147°, 
respectively. The corresponding size at 850 //m is ~ 10.0" x 7.4". The emission is centered 
on A to within 1 " in all single dish maps. The same is true for the lower S/N single-coverage 
UKT14 maps at 1.1 mm, 750, 450 and 350 /xm. In Fig. 4, we show the azimuthally averaged 
flux densities of 116293A and its surrounding envelope, and of the prestellar core I16293E. 
The intensity falls off much more rapidly for 116293A than for the prestellar core I16293E 
at radii r < 20-30 ", beyond where the slope flattens out to be similar to the prestellar core. 

The size and PA seen in our submillimeter maps agree within the errors with the C^^O 
emission mapped by Mundy ct al. (1992) with OVRO and differs in PA from the alignment 
of the two submillimeter sources A and B. In fact, neither the C^^O emission nor the 450 fim 
emission shows any clear enhancement at the B position which is the stronger continuum 
source in the wavelength range 1.3 cm to 2.7 mm. At longer wavelengths (850 and 800 //.m) 
the PA approaches that of the binary suggesting that at these wavelengths the emission 
from B is still significant. However, none of our maps lines up exactly with the binary, which 
is what one would expect if the emission origniated from a circumbinary disk. Looney et al. 
(2000) claim that their high resolution maps show that the dust emission is aligned with the 
binary. However, considering the S/N in their map, one could equally well interpret their 
result as a circumstellar disk surrounding A with a position angle similar to that seen in 
C^^O and 450 /xm, and B as a separate but partly overlapping source. Their observations 
show that about 90% of the dust emission surrounding A is spatially extended, while B is 
largely unresolved at all spatial scales they could measure. Schoier et al. (2003) observed 
the continuum at 1.37 mm with OVRO at 3" resolution and find that I16293A and I16293B 
both have a disk with diameter < 250 AU. 

We conclude that the submillimeter dust continuum and most of the molecular emission 
are centered on A and have a disk-hke morphology. The dust disk is surrounded by faint 
extended emission from the surrounding dark cloud core which falls off rather rapidly toward 
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the south and southwest. The dark cloud core is more extended toward the north and 
northeast (~ 70"), and in the west a narrow dust bridge connects it to I16293E (see Fig. 1). 

In contrast, I16293E is very extended and has a much flatter emission profile than 
I16293A (see Figs. 1, 4). The peak emission is centered on aaooo.o = 16'* 32"* 28.84^ 
(^2000.0 = —24° 28' 57.0" with a roughly triangular emission region around the peak. It has 
a ridge-likc structure at 450 fim. (see Fig. 5). The integrated fluxes (corrected for error 
lobe contribution) are 20 Jy and 133 Jy at 850 and 450 ^m, respectively, in a 50" radius 
(8000 AU). In order to perform single position photometry we also made gaussian fits to our 
submillimeter images resulting in a source size of 32" x 17" with PA = 3 ± 5° superposed 
on a more extended background. This source size was used to convert our photometry into 
integrated intensities for further analysis. 

3.2. Excitation Analysis 

3.2.1. The Protostellar Disk Around 1 16293 A 

The submillimeter maps allow us to derive a mass estimate for the dust which is more 
accurate than previous values. I16293A has been used extensively as a secondary calibrator 
for submillimeter continuum observations at JCMT (Sandell 1994), and has accurately known 
fluxes in all submillimeter bands. Although some of our UKT14 maps have relatively poor 
S/N, they still allow us to make a reliable estimate of the flux density of the compact dust 
disk, because the underlying emission from the envelope can be determined and subtracted. 
This is not possible from photometry data (Sandell 1994). In Table 1 we list the integrated 
fluxes for I16293A from two-dimensional Gaussian fits to our maps after subtraction of the 
extended envelope emission. We have not mapped the region at 1.3 mm, but instead use 
data from Mezger et al. (1992) who derived the integrated flux in a similar fashion. 

We could use the SCUBA data together with co-added IRAS data at 100 and 60 /xm 
to constrain the dust temperature, but it is clear from our submillimeter maps that the 
arcmin size of the IRAS beam will also include emission from the surrounding dark cloud, 
so to include the IRAS data we should do a two component flt and simultaneously solve 
for the dust envelope and the disk. Instead, we take a simpler approach. We solve for 
each separately and use the IRAS data to make sure that the results we get are plausible. 
Our deep 850 and 450 //m SCUBA maps predict that about half or more of the continuum 
emission seen by IRAS is likely to come from the envelope surrounding the dust disk. We 
therefore divide the IRAS 100 jim emission into halves and assume one half to originate in a 
compact < 7.4 " disk and the rest from a cooler cloud envelope around the protostar. 
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If we assume that the dust grains can be characterized by a single temperature, Tj, the 
flux density Si, at frequency u can be written as 



where is the source solid angle, B,y{Td) is the Planck function, Ti, the optical depth of the 
disk, and r^nv optical depth of the envelope. Wc assume that the envelope is optically thin 
at submillimeter wavelengths (renv <^ 1)- We write t^, as 



where Tq is the dust optical depth at frequency pq, and /3 is the dust emissivity describing 
how the dust opacity {ki, oc t,^) changes with frequency (Hildebrand 1983). We adopt the 
Hildebrand (1983) opacity (ki.2 thz = 0.1 cm^ g~^) and assume a gas-to-dust ratio of 100. 
It is now straight forward to do a least-squares fit to equation (2). Since we mapped the 
whole cloud, we know fig, and use this to constrain the fit. We omit flux densities measured 
at 3 mm with aperture synthesis telescopes, because these resolve out the extended emission 
and may also include free-free emission. 

For the submillimeter disk of I16293A wc derive a dust temperature = 40 ± IK 
(see Fig. 6), (3=1.6 and a source size of about 5". In order to make sure that these values 
are not determined by our assumed partitioning of the IRAS 100 yum flux estimate, we also 
repeated the fit omitting the 100 fira data and got similar results. Our fit to the dust disk 
is not affected by the 100 /im data because the dust emission in this source starts to become 
optically thin in the (sub) millimeter regime, rsso /xm ~ 0.42, which acts as a constraint 
for the dust temperature. This is unusual as dust emission is usually optically thin at 
(sub) millimeter wavelengths and does not provide any constraints on the dust temperature, 
unless P is known. From the fit we derive a total mass of 1.8 M© corresponding to an average 
density of n(Il2) > 10^ cm~^, and a bolometric luminosity of 16.5 L0. The uncertainty in the 
mass, due to uncertainties in the fitted Td and (3, is of the order of 0.3 M©. The continuum 
measurements of Schoier et al. (2003) show that the source size is smaller than 3". Schoier 
et al. assume optically thin emission at 1.37 mm, and a dust temperature of 40 K and derive 
a lower limit to the mass of 0.25 Mq. Our derived mass of the dust disk is in the mid-range 
of the values quoted by Mundy et al. (1986). Our fit appears to shghtly overestimate the flux 
densities at long wavelengths, which indicates that the data cannot completely be described 
by a single dust temperature. 



5, = l],5,(Td)(l-e-^'')e 



(2) 




(3) 
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3.2.2. The Prestellar Core I16293E 

Since the continuum emission traces the integrated density along the hne of sight, the 
SCUBA images can be used to derive the density distribution. We used the azimuthally 
averaged radial flux density of the continuum maps at 450 and 850 /xm to determine the radial 
density structure of I16293E. Figure 4 shows the radial emission profiles of this prestellar 
core; it is clearly seen that both emission profiles are well described by a power-law. Instead 
of the dust density we will use the H2 volume density as parameter in determining the density 
distribution. We adopt a dust-to-gas ratio of 1 : 100 and a power-law H2 density distribution 
of the form n{r) = no(r/1000 AU)"^, where uq is the H2 density at the arbitrary chosen 
radius of 1000 AU. Similarly, we adopt a dust temperature distribution following a radial 
power-law with index q and temperature Tq at r = 1000 AU. The dust emissivity is assumed 
to be of the form equation (3). 

The free parameter set of this model {uq, p, Tq, q, (3) is constrained by the total 
submillimeter flux, the radial emission proflles, and the spectral index between 450 and 850 
/im, respectively. A minimalization between the core model emission and the data yields 
best-fit parameters p = 1.0 ± 0.2 for both the 450 ixm and 850 /xm images, no = 1.6 x 10^ 
cm~^, and (3 — 2.0. The outer radius of the core was set to 8000 AU, the only value of i?out 
which allows a density profile fit with a single power- law index p. The above density structure 
yields a core mass Mcore(H2) = 4.35 M©. Prom these results we infer an isothermal dust 
temperature Td = 16 K, i.e. q — Q. We apply this temperature to the region r < 1000 AU 
while for 1000 AU < r < 8000 AU we let the temperature gradually rise to 20 K to include 
the transition of the dark core to the diffuse ISM where the temperature is determined by 
the interstellar radiation fleld. This temperature structure gives a good flt to the observed 
radial emission profiles. Note that the continuum emission is a convolution of the density 
and temperature, so the derived density and temperature structures are degenerate, e.g., a 
constant temperature Td = 25 K and uq — 1.1 x 10^ cm~^ also fits the emission profiles. 
However, the low-end T^ and high-end n fit yields a plausible physical description to this 
prestellar core and best fit to the observed HDO emision (see below) . 

The IRAS 100 /xm data were not used to further constrain the above parameter-set, 
IRAS was not sensitive to continuum emission from dust with T^ < 18 K. On the other 
hand the SCUBA emission refiects a convolution of density and temperature along the hne- 
of-sight. Thus, a cold high density region with T^<^1'&K can easily be hidden in the center 
of the core since a less dense warmer envelope component will always dominate the observed 
continuum emission. In general such a region will be hard to detect in continuum emission 
as well as in line emission. In the latter case one would need a transition of a species with 
a high critical density and low excitation temperature which is uniquely tracing the coldest 
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and densest part at the heart of the core and is not depleted. Deuterated molecules, and in 
particular the ground-state transitions of HDO, are excellently suited to trace such regions 
(see §1 and below). 

The inferred radial density and temperature structures from the dust are used to calcu- 
late the abundance profiles for the observed single dish observations of C^^O, C^^O, DCO+, 
HCO"'", H^^CO"'", and HDO. The modeling of the molecular excitation and line radiative 
transfer uses a spherically symmetric Monte Carlo method (Hogerheijde & van der Tak 
2000). The core was divided in 40 concentric shells. All shells were found to be optically 
thin at the modeled transitions of the molecules. We assume Tkin = Td throughout the core 
and adopt for each line a constant local turbulent linewidth. After convergence of the level 
populations was reached, the spectral line profile of the observed transition of each molecule 
was calculated for the appropriate beam. The calculations were done iteratively starting 
with an educated guess for the abundance and the turbulent linewidth, until a best fit of the 
modeled line profiles to the observed spectra was estabhshed. 

Our model spectra fit the observed emission profiles for all relevant molecules (Fig. 
7) and yield the following abundances for the CO-isotopomers: [C^^O]/[H2]=3 x 10~^ and 
[C^^O]/[H2]=9 X IQ-^ Adopting standard CO isotopomer ratios [CO] : [C^^O] = 500 : 1, and 
[CO] : [C^'^0] = 2500 : 1 (Wilson 1994), our derived abundances imply that CO is depleted 
by a factor of about ten with respect to the standard abundance. 

The width and strength of the observed HDO ground-state emission line is well fit by 
a [HD0]/[H2] abundance ratio of 2 x 10~^°. Note that the derived turbulent linewidth of 
the HDO line is much lower than that of the other lines. It is expected that the turbulent 
line width is roughly the same for all molecular species when observed with similar angular 
resolution if they are distributed similarly throughout the core. This is in fact the case for 
the other molecules observed with beams of about 15" (DCO^, HCO^, etc.) for which we 
obtain an average value of 0.7 ±0.1 km s^^. Although the weak HDO emission has a S/N of 
only 3, its narrowness suggests that it resides in the innermost region of the I16293E core. 
In this cold region (Tk < 16 K) the abundance of a deuterated species like HDO is expected 
to be enhanced through equation (1), while non-deuterated molecules are expected to be 
depleted (§1). A higher S/N HDO spectrum is required to confirm this. 

The best fit to the HCO+ 3-2 fine yields an abundance [HC0+]/[H2] of 1 x 10-^°. The 
H^^CO+ 3-2 emission can be reproduced for a constant abundance of 2 x 10"^^ through- 
out the core, and the DCO"*" 3-2 line fit yields an abundance of 5 x 10~^^. Note that 
[DCO+]/[Hi3CO+]=2.5, that is the DCO+ abudance is larger than that of Hi3co+. Assum- 
ing [CO]:[i3CO]=65:l (Wilson 1994), we find that HCO+ is depleted by a factor of about 
10. Table 3 summarizes the inferred abundances and local turbulent linewidths. 
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3.2.3. The Envelope Around I 16293 A 

In this section we determine a radial density and temperature structure for the YSO 
envelope using an inside-out collapse power-law density structure where the slope of the 
density and the velocity depend on the location of the collapse expansion wave. 

We assume that the dust temperature follows a radial power law. Such a distribution is 
expected for a spherical cloud with an embedded heating source in its center. Together with 
the FIR luminosity and distance of 160 pc, the temperature distribution is constrained from 
the total observed submillimeter continuum emission from the disk and the envelope (Table 
1). The inferred temperature profile ranges from — 115 K at the inner radius Tin = 100 
AU, to Td = 12 K at the outer radius rout = 7300 AU. 

We use the spherically symmetric self-similar solution of a collapsing cloud core derived 
by Shu (1977) to determine the density profile. This model is characterised by only two 
parameters: the sound speed a, and the time t since the collapse starts at the center at 
t — Q. The initial stationary density distribution is an isothermal sphere where the density is 
proportional to r~^. At a time t after the onset of collapse, the head of the collapse expansion 
wave is radius tcew — ^t- Inside this radius the velocity field varies from stationary to free 

1 _3 

fall, V oc r 2 J while the density varies as n oc r 2 . We determine the infall parameters 
independently from the molecular spectroscopic data as well as from the dust continuum 
measurements. We start with the former. 

We do not use RxA2 observations, which instrumentally broadens the lines, but only 
345 GHz window spectra to determine the infall parameters. Earlier studies often used 
emission lines from lower rotational transitions of HCO^ or H^^CO^. Our line profiles of 
HCO+ 4-3 and H^'^CO^ 4-3 are complex with excess emission at the red and blue wings of 
the doubly peaked emission profile, respectively (see below). This points to a contribution 
by the outflow(s). The red and blue wings of the 000+ 5-4 spectrum are symmetric, and 
we therefore use this line profile to constrain the parameter set a and t. These parameters 
respectively determine the width and integrated line intensity of the profile. 

A good fit to the wings of the DCO"'" 5-4 emission is found for t = 2.5 x 10^ yr, a = 0.7 
km s~^, and a constant radial abundance [DCO+]/[H2]=2 x 10~^^ (Fig. 8). This yields an 
envelope H2 mass M — 6.1 Mq within a radius of 7300 AU. Note that the absorption can 
easily be fitted by adding a cold outer shell where the temperature has dropped to < 10 
K. Our derived stationary envelope mass corresponds well with the mass derived by Schoier 
et al. (2002) of M = 5.4 Mq inside 

^out = 8000 AU, but is as much factor of four 
higher than the masses derived by Ceccarelli et al. (1996) (M = 1.88 Mq, rout = 5307 
AU) and Narayanan et al. (1998) (M = 2.3 Mq, rout = 6000 AU). In Fig. 9 we plot the 
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calculated dust and density profiles. For comparison we also plot the independently derived 
density and temperature profiles from Schoier et al. (2002) and Ceccarelli et al. (1996) on 
a log-log scale, so power-law distributions become straight lines and the outer regions of 
the envelope stand out more clearly. Our profiles agree well with those of Schoier et al. 
(2002) who derived the density and temperature profiles independently from a subset of 
our continuum measurements. Our values for the infall parameters compare well with those 
derived by Schoier et al. (a = 0.65 — 0.95 km s~^, t = (1.5-3.5) x 10^ yr) and Ceccarelli et 
al. (a = 0.5, t — 2.3 X 10^ yr). The age is more than a factor of two lower than derived by 
Narayanan et al. (1998). This is probably because our study is focused on high- J transitions, 
and therefore more sensitive to the warm inner envelope regions, while previous work dealt 
mostly with lower J transitions which are dominated by the gas in the colder outer envelope 
regions. 

We used the azimuthally averaged radial emission profile (Fig. 4) to determine indepen- 
dently the infall parameters from the SCUBA continuum maps. A best fit of the inside-out 
collapse model from Shu (1977) yields a sound speed a = 0.7 km s~^ and an age t = (0.6- 
2) X 10^ yr, close to the values derived from our fits to the line profiles. This indicates that 
the model of Shu can be used to predict the infall velocities for this class of YSOs, but 
that the age is better constrained from the line profiles. Schoier et al. (2002) also found a 
discrepancy between the infall parameters derived from dust continuum and molecular line 
emission. The infall analysis of Hogerheijde & Sandell (2000) for a sample of YSO envelopes 
indicate that this discrepancy may be a common property of Class YSOs, since a better 
correspondence was found for Class I and older sources. A possible explanation is that col- 
lapse models generally predict a much higher mass-accretion rate early in the evolutionary 
stage than the Shu model does. We therefore use the infall parameters a and t from the 
DCO^ fit to model the profiles of the observed molecular lines and to infer the abundances, 
the results are presented below. 

HDO 

The best fit to the observed HDO profile imphes a constant abundance of 3x 10~^° throughout 
the envelope. The model profile is narrower than observed. This may indicate a significant 
contribution to HDO emission from, e.g., the outflow(s) (§4.3). To reproduce the fine wings 
by infall, the sound speed must increase to a = 1 km s~^. The fit to the absorption could be 
improved by adding a cold shell with T < 10 K. But, HDO collision rates are only available 
for T > 50 K (Green 1989). Our models use the 50 K rates at lower temperatures while the 
true rates may easily be a factor of two lower at T ~ 12-1 5 K enhancing the optical depth 
of the HDO line. We therefore expect that extended collision rates would improve the fit 
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of the line profile. The continuum level results from the included FIR radiation field in the 
excitation analysis and matches the observed level very well (Fig. 8). 

We use the best HDO infall model fits {a — 1 km s~^, t — 2.5 x 10^ yr) to model 
the HDO 2 ii~2i2 and 3i2-2i2 lines. These lines have energy levels far above ground and 
are not expected to be coUisionally excited for this temperature range. However, in high 

mass YSOs these lines are found to be much stronger than expected on basis of collisional 
excitation (Jacq et al. 1990; Helmich et al. 1996b; Gensheimer et al. 1996), probably because 
the high levels are populated through excitation by the FIR radiation field. With the earlier 
mentioned dust parameters for the FIR excitation, the model yields a 2ii-2i2 peak intensity 
of = 0.006 K, and 3i2-2i2 = 0.001 K. The former value is much lower than the upper 
limit to the line strength of 0.08 K from van Dishoeck et al. (1995) for the 2ii-2i2 line. Since 
the upper-limits of van Dishoeck et al. for the 2ii-2i2 and 3i2-2i2 lines are high with respect 
to what can be achieved with present-day, more sensitive receivers, we decided to reobserve 
both lines. Deep spectra reveal no featiires at a noise rms of 0.02 K at either frequency for 
a spectral resolution of 1.25 MHz (=1.7 km s~^). 

The excitation of H2O has been calculated using the collision rates of ortho- and para-H2 
with ortho-H20 (Phillips, Maluendes & Green 1996) for kinetic temperatures between 20 and 
140 K for the lowest 5 rotational energy levels of H2O. This range traces the temperature 
structure of our model very well. We interpolated the rates for the appropriate temperature 
of each shell. The ortho/para H2 ratio in each shell is assumed to be in LTE at the appropriate 
temperature of each shell; for low temperatures, Tkin < 70 K, most of the H2 is in its para 
modification and [ortho-H2]/[para-H2]— > for lower Tkin. 

We included FIR dust excitation in our model, using the dust emissivities from Os- 
senkopf & Henning (1994) . Calculations with a large range in power- law emissivities show 
that the dust emission has not much impact on the excitation of the lowest water lines in 
the envelope of I16293A. Therefore, we did not consider the FIR radiation in our further 
calculations. The total H2O abundance is derived assuming an ortho-to-para H2O ratio of 
3 : 1. 

The abundance was varied to fit the self-reversal and to match the Iiq-Ioi peak H2O 
emission observed by SWAS. Our infall models fail to match the large width of the observed 
H2O spectrum. This is not surprising; in the large SWAS beam (FWHM^ 4 arcmin 35, 000 
AU) emission from our model cloud is dominated by the cold outer envelope. The width of 
the observed line profile clearly indicates that most of the H2O emission arises in molecular 
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outflows associated with this Class YSO which are not included in our envelope model. We 
therefore focus here on the region causing the absorption and the emission at the systemic 
velocity of the cloud, and model the H2O outflow in §4.3. Our envelope model predicts a 
narrow double-peaked profile which can match the intensity level of the observed maxima, 
but not the width (Fig. 10). Initially, the analysis was done for a uniform abundance [ortho- 
H20]/[H2]=2 X 10~^ across the envelope. This fits the maximum emission, but causes too 
deep an absorption. To reduce the absorption we applied a step function for the abundance 
profile with a drop in the H2O abundance at Tkin = 14 K, where atomic oxygen starts to freeze 
out. Such a two-step abundance profile fits the observations for [ortho-H20]/[H2]=3 x 10"^ 
for Tkin < 14 K (Fig. 10). The FIR H2O fines observed with ISO (CeccareUi et al. 2000) 
are much stronger than predicted by our model, indicating that they do not originate in the 
envelope. 

The H2D+ 

lio^lii line has been marginally detected at a 3-sigma level. A constant abundance 
of 3 X 10~^° yields a good fit to the observed maximum emission, but the fine width of 
the model spectrum is too large. We therefore tried a step function for the abundance, 
[H2D+]/[H2]=l X 10-12 foj. rj.^^^ > 20 K and [H2D+]/[H2]=2 x lO"'^ for Tkin < 20 K which 
yields a good fit to the complete fine profile (Fig. 8). 

iV2/f+ 

The N2H+ 4-3 line was observed simultaneously with the H2D+ Iiq-Ih line. A constant 
abundance with the radius of 3 x 10-^^ yields a good fit to the observed line profile (Fig. 8). 

HCO^ 

We can only put a lower limit to the HCO"*" abundance of 1 x 10"^, where the fine 4-3 
starts to become optically thick. However, this does not match the peak intensities well, 
nor does it fit the apparent extra blue and red components. The wings are well matched 
for a = 1 km An outfiow is likely contributing to the emission. The abundance for 

^kin < 10 K could be enhanced to get a deeper absorption, but this would not add much 
information since the emission and the absorption are already optically thick. Note that the 
high optical depth of the emission and absorption masks any signature of chemical changes 
at low temperatures. This transition is therefore not suited for an accurate derivation of the 
abundance structure. 
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For this HCO^-isotopomcr, we also have problems in modeling the observed double line 
structure, probably because kinematically distinct components are present at the blue side 
(Fig. 8). An abundance of 2 x 10~^^ gives a good fit to the red component. This component is 
optically thin and we attribute it to the envelope. Note that a constant abundance thoughout 
the envelope is used for this fit. The blue component is probably associated with one of the 
quadrupole outflow components (§4). 

The CO 6-5 and 7-6 lines have been modeled for a ^^CO abundance [i2CO]/[H2] = l x 10^"^. In 
Fig. 11 the model spectra are overlaid on the observations. It can be seen that the envelope 
model fails to fit the observed emission for a sound speed a = 0.7 km s~^. Even for a = 1 
km s~i, the width and maxima of the observed spectra cannot be matched. This indicates 
that, hke the H2O emission, most of the observed high-J ^^CO emission is associated with 
the outflows (§4.3). 

The C^^O 3-2 line profile can well be fitted for an abundance of 2 x 10~^ throughout the 
envelope, corresponding to a 'standard' [C0]/[Il2] abundance of about 10~^. Schoier et al. 
(2002) derive a C^^^O abundance of 6.2 x 10^^ on basis of similar observations with the 
JCMT. However their abundance yields a fit to the spectrum which is about a factor of 
two lower than our fit (see their Fig. 4), so the agreement is within a factor of 1.5. Our 
abundance-fit yields a C^^O 2-1 spectrum whose central region is quite different from the 
observations indicating a lower abundance. The observed low emission level of the C^^O 
2-1 may be caused by a low-density C^^O component with a 'standard' abundance residing 
in dark/translucent transition region between the outer envelope and the diffuse interstellar 
medium, while CO in the cold outer envelope regions may be highly depleted. Moreover, 
the modeled C^^O 2-1 and 3-2 line profiles are somewhat narrower than observed. The line 
wings of both transitions can be fitted by infall if a sound speed of a = 0.9 km s^^ is adopted, 
or there may be a contribution from the outfiow. Figure 8 plots the excitation model fit for 
the C^^O 3-2 line. We conclude by noting that the low- J C^^O transitions are not well suited 
to determine the CO depletion in the envelope since the fine profile indicates than C^^O is 
sensitive to material which may not reside in the envelope. 
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We use a constant C^^O abundance of 4 x 10~^, i.e., a 'standard' CO abundance for Tkm > 20 
K and depletion by a factor of 50 for Tkm < 20 K. This fits better than a uniform abundance 
throughout the envelope, because depletion makes the line appear broader. The line is also 
broadened by unresolved hyperfine lines, e.g., visible as a low- velocity shoulder. Still, the 
observed line is broader than our model predicts. The wings fit better for a sound speed 
a — 0.9 km s~^. There is a hint for a depression of the peak emission. Schoier et al. (2002) 
derive a C^^O abundance of about 2 x 10""^ which agrees within a factor of two with our 
results using a constant abundance value. 

The ortho-H30+ ^,^-2^ line at 396.272412 GHz has been searched for in I16293A by Phillips 
et al. (1992). We used their upper limit and our envelope excitation model to constrain 
the abundance [ortho-H30+]/[H2]< 1.5 x 10~^. This correspondes to an [H30"'']/[H2] upper 
limit of (3-5)xl0~^ if the H3O+ resides mainly in the Tkin < 50 K region, or Tkm > 50 K 
region, where the [ortho-H30+]:[para-H30+] ratio varies between 1 : 2 and 1:1, respectively. 
These upper limits are about an order of magnitude larger than derived by Phillips et al. 
(1992) using statistical equilibrium calculations for a slab with a homogeneous density and 
temperature. 

The derived abundances of all species are summarized in Table 3 together with those of 
Schoier et al. (2002). Taking into account the differences between their fits to the observed 
spectra and ours (see above), the abundances agree within a factor of 1.5. 

4. The CO Outflows and What Drives Them 

4.1. The NE-SW Flow: Powered by 116293 A 

The symmetry axis of the NE-SW outflow, as determined from our CO 2-1 map, goes 
through A at all velocities (Fig. 12). The centroid of the overlapping blue- and red-shifted 
emission is centered on I16293A, which is also evident from the velocity position plot in 
Fig. 13 (although inescapably includes some of the E-W outflow). Furthermore, aperture 
synthesis maps of dense outflow tracers, especially CS (Walker et al. 1993) and SO (Mundy 
et al. 1992) clearly show high velocity gas emanating from I16293A in the direction of the 
CO outflow, while I16293B shows no activity, nor does it seem to be associated with a 
maximum in any high density gas. We therefore attribute the NE-SW outflow to I16293A. 
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Our submillimeter maps (Figs. 1 and 4) clearly show that the dust emission associated with 
the free-free and millimeter-source I16293A is extended. This is one of the few YSOs where 
we can resolve the dust emission: it has a disk-like structure with a size of roughly 2000 AU. 
At low frequencies the continuum emission from I16293A splits up into a double source. We 
do not believe that 116293A is in itself a binary; a more plausible interpretion is that we 
are seeing highly coUimated ionized gas from the outflow, as has been seen in L1551 IRS5 
(Bieging et al. 1985; Rodriguez 1986), where the jet dominates the emission at low frequencies 
and the disk emission becomes dominant at high frequencies. These are strong arguments 
for an accretion disk that is centered on I16293A and drives a powerful outflow, which is very 
well coUimated close to the star (e.g., free- free emission), but relatively highly coUimated 
even in CO. The ratio of length to width is about 1:3 at high velocities. The outflow is at 
high inclination, since we can see overlap between the blue- and the red-shifted lobe. If we 
assume the outflow to be conical, which according to the velocity position plot (Fig. 13) 
appears to be a reasonable assumption, we can use the method outlined by Liseau & Sandell 
(1986) to estimate the inclination of the outflow, and flnd about 65° for both the blue and 
the red lobe. This is slighly higher than what Hirano et al. (2001) find on the basis of their 
SiO 2-1 map for the NE red lobe, where they derive an inclination by 40-45° from the plane 
of the sky. 

The properties of the outflow have been derived in several studies (e.g.. Walker et al. 
1988; Mizuno et al. 1990; Castets et al. 2001; Hirano et al. 2001; Garay et al. 2002; Lis 
et al. 2002, but we will redo it here; partly because we have more accurate data, but also 
to point out how uncertain these estimates are. The deep integrations of ^^CO and C^^O 
toward I16293A (Fig. 2) clearly show high velocity emission, which indicates that the ^^CO 
optical depth in the wings is substantial. At near outflow velocities (~ 1.5 km s^^ from the 
cloud systemic velocity) we flnd optical depths of 40 or more, both from ^^CO and C^^O, 
but note that here the line wings originate from both outflows and also get a signiflcant 
contribution from the underlying accretion disk. We therefore obtained additional ^^CO and 
^^CO spectra at the peak positions in the lobes of the NE-SW outflow, i.e., at offset positions 
(+80",+55") for the redshifted lobe and (-130",-40") for the blueshifted lobe. These show 
that the red-shifted outflow has substantially higher optical depth. We estimate a ^^CO 2-1 
optical depth of about unity at radial velocities as high as 11 km s"-*^ away from the systemic 
velocity of the cloud, whereas the optical depth has dropped to unity at -|-6km s"^ away 
from the cloud core for gas in the blue-shifted outflow. The ^^CO spectra are not sensitive 
enough to enable us to derive optical depths of less than about unity. 

We compute the outflow mass, momentum, momentum flux, and energy content with 
and without opacity corrections, by assuming that we can apply the same opacity correc- 
tion to all of the gas in the outflow. Since we have not directly measured the excitation 
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temperature of the outflowing gas we assume it to be 50 K, but note that it could be as 
high as 80 K (van Dishoeck et al. 1995). The values in Table 4 arc obtained assuming an 
abundance ratio 10""^ for [C0]/[H2] and summing over the map in 1 km s~^ wide velocity 
intervals starting 1 km s^-*^ from the cloud systemic velocity (assumed to be 4 km s^^). The 
apparent dynamic outflow time scales are about the same for both lobes, (3.0-3.5) x 10^ yr. 
This agrees with the dynamical timescale of (5-7) x 10^ yr derived by Hirano et al. (2001) for 
the NE red lobe. However, the blue outflow lobe appears more energetic. Normally when 
one sees an imbalance in the momentum or force of the outflow, it is an indication that one of 
the outflow lobes has penetrated the cloud surface and therefore has no gas to interact with. 
Here we see no evidence for this. A more likely explanation is that we have not properly 
corrected for the opacity and excitation differences in the two outflow lobes. The blue lobe 
is more extended and has higher apparent velocities, suggesting that the gas surrounding 
the blue lobe is less dense than the gas surrounding the red outflow lobe. This is also seen in 
our opacity estimates, which are higher for the red-shifted outflow. A similar conclusion is 
reached by Hirano et al. (2001) and Garay et al. (2002) on basis of SiO and CH3OH maps, 
respectively. 



4.2. The E-W Outflow: A Fossil Flow Driven by I16293B 

The symmetry axis of the E-W outflow passes north of I16293A, but not as clearly 
deflned as the NE-SW outflow. The red-shifted outflow lobe is much more compact than 
the blue-shifted lobe, and appears to bend more towards the north. However, the peak 
intensity of the outflowing red-shifted gas is rather well ahgned with the large blue-shifted 
outflow lobe with an average position angle for the flow of PA ~ 96 ± 3°. The E-W outflow 
appears to intersect with the continuum disk about 5-10" north of I16293A. Since there is no 
source near this position other than I16293B, we associate the E-W outflow with I16293B. 
A similar conclusion was also drawn by Walker et al. (1993), although they stress that at 
present there is no evidence for outflow activity in the immediate vicinity of B. Garay et al. 
(2002) flnd that the CH3OH wing emission has a narrow and nearly constant velocity width 
along the E-W lobes, which they attribute to turbulent entrainment. This is consistent with 
our flnding that the E-W outflow is a fossilized flow. 

It has been suggested that the NE red outflow lobe and the eastern blue outflow lobe 
could be part of a wide angle flow similar to the L 723 outflow (Walker et al. 1988) or 
that the two outflows arc part of a precessing or episodic outflow driven by the same star 
(Walker et al. 1988, 1993; Mizuno et al. 1990). However, our CO map clearly shows that 
these are two separate outflows. Since the symmetry axis of the E-W outflow does not 
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intersect with I16293A, it is rather unhkely that this outflow would have been caused by 
an earher outflow episode from I16293A. Hirano et al. (2001) and Garay et al. (2002) also 
conclude that two independent bipolar outflows are present on the basis of their SiO and 
CH3OH maps. Mizuno et al. (1990) also discuss a third outflow, because in their low spatial 
resolution Nagoya map, the E-W outflow extends more than 10' from IRAS 16293—2422 and 
their high resolution Nobeyama map does not show this extended outflow component. Our 
map does not extend as far as the Nagoya map, and our spatial resolution is poorer than that 
of the Nobeyama CO 1-0 map. However, the CO map we present here has the advantage of 
being fully sampled, unhke the SiO 2-1 emission mapped by Hirano et al. (2001), and it has 
a much higher resolution (20") than the large-scale SiO 2-1 map of Garay et al. (2002) (57"). 
We therefore better resolve the two outflows. One can clearly see that the outflow continues 
to the east outside the area we have mapped, but with rather low velocities. Therefore the 
extended blue low-velocity outflow appears simply a continuation of the E-W flow. 

We also clearly see the interaction of this outflow with the I16293E core, which was 
already found by Mizuno et al. (1990). In our channel maps, which are overlayed on the 
800 ;um-continuum emission (Fig. 12), this interaction is clearly visible. It appears rather 
gentle: the outflow is enhanced at low velocities and seems to stream around the core and 
forms almost a cavity immediately behind it. That it streams around the core region is 
supported by the fact that we see faint red-shifted ^^CO towards the dense core region, but 
the red-shifted emission is much stronger in front of the core, although one can still see it 
faintly on the backside as well (Fig. 12). Figure 5 shows a blow-up of the 450 iim image of 
I16293E. It is clear that the outflow has compressed the core and caused the high density 
N-S ridge that we see in our continuum maps and prominently at near cloud velocities in 
the channel-maps of Fig. 12. Behind the core it appears that the outflow gets recollimated. 
There is no indication of interaction between I16293E and the red-shifted NE outflow lobe 
from our CO channel maps, nor from the SiO channel maps (Hirano et al. 2001; Castcts et al. 
2001; Garay et al. 2002). Castets et al. (2001) attribute the blue and red stream of the E-W 
outflow around the I16293E core to the red and blue lobe of an outflow powered by I16293E 
causing two shocked regions (HEl and El in their nomenclature) visible in their H2CO maps. 
They support their hypothesis by means of HCO+ 1-0 spectra taken along a shce which goes 
through I16293E as well as through the peaks in the H2CO emission, since these spectra 
reveal a blue wing at one side of I16293E and a red wing on the counter side. Although 
these wings are real, the results are different when put in context of the entire cloud. Our 
large maps show what small maps, which do not fully cover the complete interaction region 
of the extended outflows with the LDN 1689N cloud, cannot: local small scale flows in the 
vicinity of I16293E are part of large scale motions. From our CO map it is clear that the 
apparent SW-NE outflow and the alignment with I16293E are due to projection effects of 
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pieces of the E-W outflow from I16293B which streams around the dense prestellar core 
I16293E. This view is supported by the study of Hirano et al. (2001) who find that the E-W 
flow has fan-shaped lobes, each of which containing blue- and red-shifted components. This 
would imply that El is red-shifted gas in the far side of the eastern blue lobe. The red-shifted 
shocked region El is thus the interaction region of the E-W outflow, while the blue-shifted 
H2CO peaks HEl and HE2 from Castets et al. (2001) are the interaction of the near side of 
the El-W blue lobe with dense gas of LDN 1689N located at the southwest and southeast 
side of the dense core I16293E. These interaction regions could also be the edges of I16293E. 



4.3. The High-J CO, HDO, and H2O Emission: Outflow Modeling 

Section 3.2.3 found that the envelope model of I16293A can reproduce the line profiles of 
most observed molecular transitions, but not of CO 6-5 and 7-6, HDO, and H2O. Although 
the absorption at the cloud systemic velocity can be fit with an envelope model, the broad 
emission cannot be reproduced even for sound speeds as high as a = 1 km s~^. Here we 
investigate the origin of the broad high- J CO, HDO, and H2O wings by means of an outfiow 
model. 

We have scaled our CO 7-6 emission to the level of the H2O emission (Fig. 14). It is 
striking that the two spectra have such a similar shape despite the fact that the beam sizes 
differ by a factor of about 30. The match of the absorption widths is also very good, and 
indicates that we may see here the same region with JCMT and SWAS in absorption. In 
addition, the line wings are similar, although at the blue side there is an extra CO emission 
component while at the red side there is an extra H2O emission component present. But 
the extent of the wings is similar for both species. Does this mean that the wings are from 
a very small region covered by both beams, or that the wing component is rather uniformly 
distributed throughout the cloud? 

The SWAS data were fitted with a Monte Carlo radiative transfer code and a physical 
model of a spherically symmetric cloud subdivided into shells of constant physical conditions. 
The use of the code and its application to water lines is described in Ashby et al. (2000), 
which also analyzes the variety of water line profiles and establishes broad diagnostics of 
outfiow, infall, and turbulence apparent in SWAS data. The fine shape lead us to explore 
models which combine outfiow and turbulence; we find that the line strengths depend on the 
velocity field as well as on the temperature and abundance in a complex relation. Figure 10 
(middle panel) shows the predicted ortho-H20 spectra from the outfiow model with an ortho- 
H2O abundance of 1 x 10^^ throughout, a turbulent velocity of 2.3 km s^^, and maximum 
outfiow velocity V = 5.8 km s~^. The temperature ranges from 200 K to 15.5 K. The outfiow 
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model provides an acceptable fit to the observed extended wing emission, but cannot explain 
the central part of the observed profile. The outflow model results are complementary to 
the infall modeling of the ortho-H20 line (Fig. 10, top and middle panel), which fits only 
the central region of the observed spectrum. To illustrate this qualitatively, we plotted in 
Fig. 10 (bottom pannel) the sum of the outflow and envelope spectra of ortho-H20 which 
indeed yields a rather good description of the entire SWAS spectrum. Note that the red 
wing of the H2O emission is not fitted by our model, this may be due to multiple outflows 
(§4.1 and §4.2) which we did not model. For a quantitative comparison we do not only need 
a mutiple onflow model but also have to incorporate the outflow and infall models in one 
single excitation model since they overlap in extent and velocity space. This could mean 
that part of the blue infall gets absorbed by the blue outflow component and vice versa. An 
avenue for further study would be a combined infall and outflow model. 

The same outflow models were used to model the CO 6-5 and 7-6 lines in the JCMT 
beam using an abundance [CO]/[H2]=10~'^ (Fig. 11). The wing emission predicted by the 
outflow model poorly fits the observed emission which indicates the presence of multiple 
outflows, analogous to the H2O situation. 

Finally, we calculated the HDO emission in the outflow in a 11" beam. We started with 
a uniform [HD0]/[H2] abundance of 3 x 10^^'^, similar to the envelope value. This yields a 
spectrum which matches the maximum intensity of the observed emission, but not the wings 
(Fig. 15). This may be due to the presence of multiple outflows, similar to the H2O case. 
The outer wings can be fitted for an abundance of 3 x 10~^ and results in a self-reversed 
profile, but its peak intensity is too high by a factor of 2-4. This is not surprising since we 
do not include an infalling envelope. A better fit may be obtained by a combined (multiple) 
outflow and infall model as discussed above for H2O. We conclude that our shock model 
can explain the wings of the observed HDO emission and that the HDO abundance in the 
outflow may be somewhere in the range from 3 x 10~^° to 3 x 10~^. 



5. The [HD0]/[H20] Abundance Ratio Structure Toward I16293A 

5.1. Derived [HD0]/[H20] Profile From Envelope Excitation Analysis 

It is clear that great caution should be taken when comparing the H2O abundance 
derived from the SWAS measurements and the HDO abundance derived from the JCMT 
observations. Nevertheless we arc able to use the same models since we convolve the emerging 
emission with the appropriate beam. Because of the large bcamsize, SWAS is not sensitive to 
the H2O emission from the warm inner envelope. To test this we enhanced the abundance of 
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H2O for Tkin > 90 K to the value of 10~^. This was not found to have any significant impact 
on the emission in the SWAS beam. Clearly SWAS is mostly sensitive to H2O emission from 
extended colder gas in the outflow, and to a lesser extent to gas in the outer envelope region 
(§4.3). 

The shocks associated with the outflows are able to release H2O from the grains and 
effectively convert oxygen into H2O through neutral-neutral reactions, but the velocities are 
so small that H2O formation at an 1 x 10"*^ abundance level is inhibited. Indeed, ISO LWS 
observations yield much lower H2O abundances. The ISO observations of H2O are all from 
higher level transitions E/k > 100 K, and thus trace the hot/warm disk/core region of 
small extent. Ceccarelli et al. (2000) derive an [H20]/[H2] abundance in the 'outer envelope' 
(r > 150 AU) of about 2.5 x 10"^ which is similar to our envelope abundance (2 x 10~^), 
while for the inner region (r < 150 AU) they derive an enhanced abundance of 1.5 x 10"^ 
which is attributed to evaporation of water ice mantles of grains at Tkm — 100 K. The latter 
value is higher than our inferred outflow abundance of about 10~^ (§4.3). 

If we combine the H2O abundances of Ceccarelli et al. (2000) with our derived HDO 
abundances we find that the [HD0]/[H20] abundance ratio varies between 2 x 10~^ in the 
warm envelope, r < 150 AU, 10~^ in the envelope with r > 150 AU and Tkm > 14 K, and 

8 x 10~^ in the outer envelope where Tkm < 14 K. For a detailed comparison between the 
HDO and H2O abundances throughout the envelope we consider in the next section the 
temperature effects of ion-molecule chemistry and the effects of depletion of molecules and 
atoms which play a vital role in the formation of HDO and H2O. 



5.2. A Simple Chemical Model 

In ion-molecule chemistry (e.g., Pineau des Forets et al. 1989; Millar et al. 1991) HDO 
is thought to be formed through the dissociative recombination 

H2DO+ + e- ^ HDO + H, (4) 

with a statistical branching ratio of | compared to the corresponding recombination of HaO"*". 
In steady state the gas-phase density of H2DO"'" can be written as 

m T^n+^ _ + n(H20)A;2)n(H2D+) + n(H20)n(DCO+)A;3 

where ki are the relevant reaction rate coefficients. We adopt equal total rate coefficients 
for the reactions involving the deuterated and hydrogenated species and assume statistical 
branching ratios. HaO"^ is formed through similar reactions with hydrogenated species. We 
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will neglect the formation of H3O+ from channels leaving the H2DO+ reaction sequence and 
vice versa, and neglect the formation of H3O"'" and H2DO''' from reactions of or H2D''' 
with HCO"*", and DCO"*", respectively. The [HD0]/[H20] abundance can then be written 

x(HDO) _ 2x(H2DO+) 



x(H20) 3 x(H30+) 

1 a:(H2D+) 
18 xiRt) ' 



(7) 



The heart of the fractionation of dcutcratcd molecules is the molecular-ion H2D+ which 
is formed through reaction (1). In equilibrium the formation and destruction rates of H2D"'" 
are balanced which leads to (Stark et al. 1999) 

3:(H2D+) ^ x(HD)fcj + 3:(D)fcD 

where ke is the dissociative recombination rate coefficient of H2D+, /cd is the rate coefficient 
for the formation of H2D+ via the reaction H3 + D, x(X) is the fractional abundance of 
species X=CO, O, H2O, etc., and kf and kr are the forward and backward rate coefficients 

of reaction (1), see Stark ct al. (1999) for references of the adopted rate coefficients. The 
density depends on the effective cosmic ray ionization rate (eS and can be written as (Lepp 
et al. 1987) 

- E&y 

We consider only the destruction of H2D+ through reactions with CO, assume a;(HD) = 
10a;(D) = 2.8 x 10^^, and (cs = 5 x 10^^^ s^^. Together with the balance for formation 
and destruction (Eq.[9]) we derive an expression for the H2D+ abundance with the CO 
abundance as the only free parameter. This abundance structure was subsequently used as 
input in the excitation calculations to compute the spectrum of the ortho-H2D+ ground- 
state transition. A good fit to the observed H2D+ spectrum is established for a two step CO 
abundance profile with [CO]/[H2]=1.5 x 10"^ for Tkin > 22 K and 3 x 10"^ for Tkm < 22 K 
(Fig. 16) applied to the infall envelope model. This indicates a CO-depletion by a factor 
of 50 at temperatures below 22 K, where CO is frozen out. This value equals the depletion 
derived from our C'^^O 2-1 excitation analysis, but is a factor of 50 larger than derived from 
our C^^O 2-1 and 3-2 observations which indicated no depletion (§3.2.3). However, the 
spectra of the latter transitions are biased toward the warm inner regions and the outflow. 

Now we can use the [H2D+]/[H3 ]-model fit, together with the H2O abundances from 
the excitation analysis and the expression of the [HD0]/[H20] ratio (Eq. [7]) to determine 
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an HDO model abundance profile and to calculate a model spectrum analogous to H2D+. 
The model spectra are shown in Fig. 16. When we adjust [H20]/[H2]=1.25 x 10"''' for 
Tkin > 14 K and 2 x 10~^ for Tkin < 14 K we get a fit to the observed emission peaks. 
Note that these abundances agree within the uncertainties with the derived H2O excitation 
results. The HDO chemistry model spectrum (Fig. 16) resembles the HDO excitation infall 
spectrum (Fig. 8) very well. The fact that the absorption does not go down to zero is a 
result of the uncertainty in the HDO coUision rates for T^in < 50 K. Note that the HDO 
spectrum is calculated self-consistently from the theoretical [HD0]/[H20] expression and 
the H2O 2- value abundance profile from the excitation analysis, i.e., it does not provide an 
independently determined HDO abundance structure. 



5.2.1. Observational Results Versus Chemistry Models 

The upper- limits on the ortho-HsO"*" abundances derived from the observations by 
Phillips et al. (1992) in §3.2.3 can be used to determine an upper limit to the H2O abundance. 
Prom the formation and destruction balance of H2O, we derive [H20]/[H30+]~ 1000-1500 
where we used the experimentally determined branching ratio of 0.5 from Jensen et al. (2000) 
to form H2O in the dissociative recombination of H3O+. Furthermore we only considered 
the formation of H2O via the Hj]" channel and assume that 10% of all ions cause destructive 
dissociation of H3O+, similar to Phillips et al. (1992). The upper- and lower-limits corre- 
spond to the Tkin < 50 K and Tkm = 80 K ranges, respectively. The observations then yield 
[H20]/[H2]< 5 X 10~^ for both temperature ranges and are consistent with the abundances 
derived by Ceccarelli et al. (2000) from ISO observations. 

The derived H2O abundance in the outflow is 1.3 x 10~^ and agrees with the dynam- 
ical/chemical model of Bergin etal. (1999) for a pre-shock chemistry which evolves until 
t — 10^ yr. The same model predicts significant deuterium fractionation in the pre-shock 
gas after 10^ yr: [HD0]/[H20]~ 10"^, in agreement with our derived value for the envelope 
region where r > 150 AU and Tkin > 14 K. For the warm inner-envelope region r < 150 AU, 
we derive [HDO]/[H20]=2 x 10~^. Our gas-phase chemistry model predicts values about two 
orders of magnitude lower for this region. We attribute the observed enhancement to the 
shock induced release of solid water (HDO, H2O) from the grains. 



-29- 



5.2.2. [HD0]/[H20] in the Cycle of Star- and Planet Formation 

Our derived [HD0]/[H20] ratio in the cold (Tkin < 14 K) envelope regions is 0.08. This 
ratio is however an upper limit due to the fact that the HDO collision rates are only available 
for Tkin > 50 K. For lower temperatures we used the 50 K values, but the rates are expected 
to decrease for lower temperatures, so that the true HDO abundance may be lower than 
we derived in §4. Since oxygen starts to freeze out at Tkin = 14 K we may also expect a 
drop in HDO which is similar to H2O (§5.2. Thus, the [HD0]/[H20] ratio may be in the 
range lO'^-S x 10"^ at T^n < 14 K. The [HD0]/[H20] chemistry model (Eq. [7]) reaches a 
maximum value of ~ 0.06 when [H2D+]/[H-|"]— > 1. Thus, the observed [HD0]/[H20] ratios 
can be understood from pure gas-phase deuterium chemistry at low temperatures. 

The [HD0]/[H20] abundance ratio in the warm inner-envelope (2 x 10~^) of I16293A 
is close to the value that has been measured in the solar system: 3.0 x 10~^ in the comets 
Hale-Bopp, Halley, & Hyakutake (Eberhardt et al. 1995; Bockelee-Morvan et al. 1998; Meier 
et al. 1998). and 1.5 x 10"^ in the Earth's oceans (Lccuyer et al. 1998). These values are all 
so close as to suggest that the [HD0]/[H20] abundance ratio may be a constant in the cycle 
of star- and planet formation. The fundamental question is then: where and when is the 
[HD0]/[H20] ratio determined? There are currently two favoured scenarios (e.g., Aikawa 
& Herbst 1999; Aikawa et al. 2001): 1) deuteration from ion-molecule chemistry in the 
protoplanetary disk; 2) deuteration determined by gas phase and/or solid state chemistry in 
the dense interstellar medium, namely in the parent molecular cloud before star formation 
begins. Aikawa et al. (2002) modeled the [HD0]/[H20] ratio in accretion disks models and 
found values of the order of 10^^ for R > 100 AU in the cold midplane where Tkin — 15 K. This 
may indicate a large depletion of H2O while HDO is still in the gas phase. Unfortunately 
the JCMT and SWAS obervations cannot be used to determine the (deuterated) water 
abundances in the circumstellar disk. For the moment we can only make a comparison for 
the HDO abundance in the warm envelope. Our derived [HD0]/[H2] abundance ratio of 
3 X 10~^° is comparable within a factor of a few to the protoplanetary disk model of Aikawa 
et al. (2002) with mass accretion rate 10~^ Mq yr~^, and supports the hypothesis that the 
HDO abundance in the disk and warm envelope are equal. 

6. The nature of I16293A, B and I16293E 
6.1. I16293A 

I16293A is a Class source, i.e., is believed to be a protostar according to the classi- 
fication proposed by Andre et al. (1993). It is invisible in the near- and the mid-infrared. 
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it is associated with free-free emission so it has already formed a stellar-hke core, it drives 
a powerful molecular outflow, and it is surrounded by an envelope of dust and gas. Even 
though the temperature that we derive for the dust disk surrounding I16293A is higher than 
the canonical value of 30 K proposed by Andre et al. (1993), the ratio of submillimeter to 
bolometric luminosity is typical for Class sources, namely > 5 x 10"^. For 116293A we 
derive Lsubmm/Lboi > 2 x 10~^. Note that this is really a lower limit, since our estimate for 
the bolometric luminosity also includes the luminosity of I16293B, which is hkely to be a 
few Lq, while its contribution to the submillimeter luminosity is negligible. Bontemps et 
al. (1996), in their study of outflow properties of a sample of Class I and Class sources, 
find that Class sources in general have a much higher outflow efficiency (defined as the 
ratio of momentum flux to radiative force) than Class I sources, with an average ratio for 
Class sources of about 1000. For the 116293A outflow wc flnd FcoZ-^rad ~ 730, which is 
close to the average value found by Bontemps et al. (1996). They also included 116293A in 
their sample, but the outflow momentum they derive is significantly larger because they also 
include the E-W outflow, which is driven by 116293B. 



6.2. I16293B 

Although it is most likely that 116293A and 116293B originate from the same collapsing 
core, 116293B remains much more of an enigma. Our observations support the suggestion 
of Walker et al. (1993) that I16293B drives the large E-W outflow. Our observations also 
show that the submillimeter dust emission peaks on I16293A and more or less follows the 
C^^O emission mapped by Mundy et al. (1990). Single dish observations have insufficient 
resolution to show any contribution from I16293B at all, even though we know from VLA 
and 3 mm aperture synthesis studies (e.g., Mundy et al. 1992) that the spectral index of 
I16293B is much steeper than that of I16293A and indicative of dust emission. 

What is I16293B? Based on the scarce evidence that we have, I16293B looks hke a 
young T Tauri star, which has driven an extremely powerful large outfiow which extends 
~ 0.45 pc on the blue-shifted side. Such outfiows from low-luminosity young stars are not 
uncommon; they are seen both as large scale optical jets (Bally & Dcvinc 1997) and large 
scale CO outflows (Padman et al. 1997). 116293B is a low-luminosity star, < 5 Lq, associated 
with compact dust emission; the true luminosity may be found by means of high resolution 
millimeter-wave observations. 
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6.3. I16293E 

We derive a bolometric luminosity Lboi = 2.3 L0 for the prestellar core I16293E, and 
find that its FIR continuum spectrum resembles that of the Class source I16293A (Fig. 4). 
Castets et al. (2001) quote a dust temperature = 24 K for this source, higher than our 
value of Td — 16 K since we correct for the overlapping envelope of I16293A at the position 
of I16293E. This correction is important since both cores are connected with a bridge-hke 
structure (Figs. 1, 5). 

The temperature and density structure of I16293E resembles the conditions in the outer 
envelope of I16293A. It is very interesting that the HDO abundance of the warm envelope of 
the YSO, and the core of I16293E are similar. This suggests that some of the gas species in 
the warm envelope may still have the abundance that they had in the prestellar stage, the 
difference being that the gas in the warm envelope is excited and thus easier to detect. A 
comparison between I16293E and I16293A could therefore be very important to understand 
the formation of a first hydrostatic core, a so-called Class —I protostar (Boss & Yorke 1995). 
Such protostars lie between the prestellar clouds and the Class protostars, and occur when 
a collapsing cloud first becomes optically thick, heats up and reaches a quasi continuum 
stage. This stage is short hved, about 2 x 10^ yrs, after which a second collapse takes place 
and the final protostar is formed. Class —I protostars are characterized by a central core of 
about 200 K, embedded in a cold core of about 10 K, although the models of Boss & Yorke 
mostly reached a central source with a much lower temperature of 20 K. 

Is a first protostar present in I16293E? Most of the observed molecular hne profiles from 
I16293E show two peaks with a stronger redshifted peak, characteristic of expansion motions. 
In §4.2 we argue that this may be due to projection effects of pieces of the E-W onflow from 
I16293B which streams around I16293E. We also do not find an outflow associated with 
I16293E in our CO maps (§4.2), so it should be younger than a Class object. A central 
Class —I protostar with Tkin = 20-200 K may be present but very hard to detect in deep 
mid- or far infrared continuum imaging since the spectral signature will be only a weak 
shoulder at the Wien side of the energy spectrum of the core. Such a first protostar will be 
very compact and have little mass. Mi < 0.04 Mq, which can easily be hidden in the dense 
cold I16293E core. Detection with SCUBA seems impossible because of the beam size and 
because the outer envelope of such a cold core will be warmer due to the interstellar radiation 
field and will dominate the continuum emission. Our HDO emission, which we think comes 
from the heart of the core, has a narrow linewidth {AV = 0.6 km s^^) and om models 
yield a kinetic temperature of about 16 K. This indicates that a first protostar must still 
be rather cold at the low end of the temperature interval of Boss & Yorke (1995). Sensitive 
high resolution submillimeter interferometry of the continuum with, e.g., the Atacama Large 
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Millimeter Array may reveal a cold accretion disk. 

6.4. The Extreme Deuteration of LDN 1689N 

We think that LDN 1689N is a remarkable low-mass star-forming region, but not because 
it 'harbors a young protostar, two molecular outflows and a cold nearby molecular core' (Lis 
et al. 2002). Multiple outflows are common around YSOs and have been reported, e.g., 
in L1448 (Barsony et al. 1998; Wolf-Chase et al. 2000), HH 24 (Eisloffel & Mundt. 1997), 
and NGC 1333 (Knee & Sandell 2000). Deep studies of embedded YSOs in low-mass cores 
revealed the presence of apparently starless cold companion cores with masses of about a few 
Mq (Motte et al. 1998; Hogcrhcijdc k Sandell 2000; Shirley et al. 2000). Strong deuteration 
is expected in cold dense prcstcllar cores which arc on the verge of forming a first hydrostatic 
core because of the extreme depletion of molecules like CO which may result in an abundace 
ratio [H2D+]/[H3]> 1 and a boost of the [D]/[H] ratio in molecules (§1). Indeed, several 
prestellar clouds have now been found to have a large depletion of CO (e.g., L1544, Caselli 
et al. 2002; L68, Bergin et al. 2002) and a large H2D+ abundance (CaseUi et al. 2003). But, 
what really is remarkable in LDN 1689N is that it can be considered as a 'deuteration candle' 
with the strongest deuteration known so-far in low-mass YSOs. What causes the extreme 
deuteration in the envelope of I16923A and the core of I16293E? The answer may lie in 
the fact that I16293AB is a binary system where the YSOs and their outfiows differ in age. 
I16293B is probably the oldest of the two YSOs and the source of the fossilized E-W outfiow 
(§4.2). This outflow is responsible for the dynamical interaction of a part of the eastern blue 
lobe with a dense dark core, the precursor of I16293E. This interaction caused a shock with 
velocity Vg — 8-10 km s~^ (Lis et al. 2002) as measured at the deuterium peak D north 
of I16293E and possibly higher at the first place of impact west of the core. This could 
have released deuterated species which were condensed on the dust grains, and compressed 
the gas, which then cools efficiently to low temperature (Lis et al. 2002). Since the E-W 
outfiow causes a less energetic shock than the NE-SW outfiow (Lis et al. 2002; Garay et al. 
2002) and SiO emission is found only associated with the NE-SW outfiow (e.g., Garay et al. 
2002), we expect only deuterated species to desorb since these were frozen out relatively late 
and thus are located near the outer surface of the grain mantles. These slow shocks of the 
outflow leave the grains mostly intact, with a post shock composition like that of a grain in a 
^kin = 20 K gas. The release of deuterated species, combined with the gas-phase deuterium 
chemistry at low temperatures, boosts the deuterium enrichment to extremely high values 
just after the passage of a shock; these deuterated species freeze out again in about 10^ yrs. 

Compression from this shock resulted in an elongated N-S structure, visible in the 450 



-33- 



/im image (Fig. 5). At the north side this filamentary structure bends to the northeast 
where the gentle flow passes around the core; here the deuterium peak is located. This shock 
would cause a blue-shift of the systemic velocity in I16293E with respect to the parental 
LDN 1689N cloud. Indeed, the HDO spectrum toward I16293E is blue-shifted by about 0.3 
km s~^ with respect to the systemic velocity of the LDN 1689N envelope toward 116293A 
(Table 1). Lis et al. (2002) also found such a blue-shift of deuterated molecules toward the 
deuterium peak D which is located about 10 arcsec northwest of the continuum peak of 
I16293E. We believe that at this position the shock has the greatest chemical impact. 

This strong deuterated core may have been shocked again at a later time by the NE-SW 
outflow from I16293A causing a velocity red shift of 0.3 km s~^. However, only the non- 
deuterated species, which predominently reside in the outer core region, are red-shifted with 
respect to their deuterated forms found mainly in the heart of the core. This is exactly the 
velocity shift which Hirano et al. (2001) would expect from an interaction of the northeast 
red lobe. However, they did not detect this, since they measured only the H^^CO"*" 1-0 line 
at the position of I16293E which has a velocity equal to the cloud systemic velocity. Under 
our hypothesis this has to be measured with respect to the reference frame of the first shock. 

This second shock boosted the gas-phase deuterium chemistry to its maximum and 

the temperature is further lowered due to further compression. We think that it is this 
firecracker effect which leads to the extreme dcutcration of LDN 1689N in general, and 
116293E in particular. Multiple shocks may be common in low-mass YSOs; Larsson et al. 
(2002) suggest that Serpens SMMl is infiuenced by multiple low-velocity shocks driven by 
the outfiow. 



7. Discussion and Conclusions 

We have analyzed our submillimctcr dust and continuum data in terms of an cxcitatiom 
analysis using a model for the temperature and density structure of the disk, envelope, and 
outfiow of 116293A, B, and E. The spectral lines of H2O, HDO, and H2D"'" were subsequently 
compared with predictions from a simple chemical network. 

Figure 17 shows a schematic overview of the IRAS 16293-2422 region showing its distinct 
physical regions (disk, warm dense core, warm envelope, cold envelope, outfiow regions, etc.). 
For each region we indicate the derived physical parameters, and their molecular abundances. 
For the submillimeter disk of I16293A we derive a dust temperture = 40 K, a source size 
of about 800 AU, a total mass of 1.8 M0 and a luminosity of 16.5 L© which corresponds to a 
gas density n(H2) > 10^ cm~^. The envelope of I16293A can be well described by a collapsing 
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envelope model with parameters age t — (0.6-2.5) x 10^ yr and sound speed a — 0.7 km 
s~^. In fact the center of the expansion wave, rcEW ~ oi, is 890-3700 AU from the centre 
which corresponds to 6-23" at the distance of 160 pc, comparable to the beam size of the 
SCUBA continuum and submillimeter line emission of the species we used in the analysis. 
The outer radius of the envelope is about 7300 AU and the temparature ranges from 115 K 
in the inner region to 12 K in the outer region. The envelope mass is Menv(H2) =6.1 Mq. 
We find that for a YSO such as I16293A the line emission is beter suited than the continuum 
emission to constrain the infall parameters. The excitation analysis of the envelope yields a 
constant HDO abundance of 3 x lO"^*^ throuhout the envelope, and an ortho-H20 abundance 
of 2 X 10~^ for the regions with Tkin > 14 K, while the abundance is lowered by a factor of 
150 in the colder regions. Prom the C^''0 analysis we find that CO is depleted by a facot of 
50 for the envelope regions where T^in < 20 K. The line wings of the CO, HDO, and H2O 
emission can qualitatively be explained by a single outflow model with a turbulent velocity 
of 2.3 km s~^ and where the maximum outflow velocity is about 6 km s^^. The excitation 
analysis of the outflow model indicates that the ortho-H20 abundance in the outflow is about 
1 X 10~^ and that the HDO abundance in the outflow is in the range of 3 x 10""^ to 3 x lO"^'^. 
A more quantitative comparison requires a multiple outfiow model. We find that only two 
outfiows are present, a NE-SW outfiow with a dynamical time scale of about (3.0-3.5) x 10^ 
yr which is powered by I16293A a Class source, and an E-W outflow which is a fossil 
flow with a dynamical time scale of (6.4-7.6) x 10^ yr driven by I16293B which looks like a 
young low-luminosity T Tauri star. The prestellar object 116293E is well described with an 
isothermal core of T^in = 16 K, a core radius of 8000 AU and a power law density structure 
of the form n(r) = no(r/1000 AU)~p where Uq = n{R = 1000 AU) = 1.6 x cm'^. The 
total mass mass of this core is Mcore(H2) = 4.35 Mq. Our excitation analysis of the core 
indicate a depletion of CO and HCO+ by a factor of ten relative to the 'standard abundances' 
and that the DCO+ abundance is larger than that of H^^CO"^. Our tentative detection of 
HDO indicates the presence of a cold condensation in the heart of this core where the HDO 
abundance is about 2 x 10~^°. Our detailed CO mapping shows that I16293E has no outflow, 
so its evolutionary stage is younger than Class 0. We argue that this core may hide a flrst 
hydrostatic core, a so called Class —I object which is deeply embedded in a largely unaffected 
cold core. 

The reason for the extreme deuteration in the L16293N cloud and in particular in 
116293A, and 116293E may lie in the fact that 116293AB is a binary system where the YSOs 
and the outflows differ in age causing a flrccracker efffect of the deuteration. In this scenario 
the older E~W onflow from I16293B is responsible for the formation of I16293E out of a 
dense dark core through a slow shock where only deuterated species were desorbed from the 
grains which, in combination with the low temperatures, amplifies the deuterium enrichment 
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in the gas phase before a possible re-absorption on grains. This strong deuterated core may 
have been shocked again an a later time by the younger NW-SW outflow from I16293A 
further boosting the gas-phase deuterium chemistry to extreme values. 

It is remarkable that the abundances of many molecules in the prestellar core I16293E 
and in the warm part of the envelope of I16293A are comparable. In particular, the HDO 
abundances are similar which indicates that the warm gas in I16293A may have gone through 
the same cold pre-collapse phase that I16293E is currently in. A simple chemical gas-phase 
modeling of the deuterium chemistry requires that CO is depleted by a facor of 50 for T^in < 
22 K in order to reproduce the observed H2D+ emission from our chemical network applied 
to the envelope (T, n) structure, and that an H2O abundance of 3 x 10"'^ for Tkin > 14 K and 
4 X 10~^ for Tkin < 14 K is required to reproduce the observed spectrum from the chemical 
model. The [HD0]/[H20] abundance in the warm inner-envelope of I16293A (2 x 10~^) is 
close to the value that has been measured in the solar system and supports this. Our chemical 
modeling indicates that the [HD0]/[H20] ratio is determined in the cold gas phase prior to 
star formation, and that this ratio is conserved in the early stages of low-mass star through 
freeze out on dust grains formation and that the molecules are released unprocessed to the 
gas phase after the formation of a low-mass YSO. If this warm deuterated gas is subsequently 
locked-up in a proto-planetary disk before the high temperature gas chemistry causes a de- 
deuteration, the [D]/[H] ratio of HDO in low- mass pre-solar systems would remain conserved. 
To investigate this scenario in detail, sensitive follow-up observations of the HDO and H2O 
ground-state transitions in a large sample of low-mass YSOs are needed at high spatial and 
spectral resolution (i.e., interferometry from high altitudes [HDO] and from space [H2O]) 
and allow a study of the [HD0]/[H20] abundance evolution from prestellar cores to planets. 
The H2^0 ground-state transition may be more suitable than the ortho-H20 transition since 
the ground-state of the latter molecule is at AE/k = 34.3 K above ground, while that of 
H2*^0 is at K, and the main isotopic form of water is expected to be mostly optically thick 
in its ortho and para ground-state lines. 

The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre, on 
behalf of the Particle Physics and Astronomy Research Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the National Research Council of 
Canada. It is a pleasure to thank the JCMT staff and the MPIfR Division for Submillimeter 
Technology for their outstanding support. The Canadian Astronomy Data Center is operated 
by the Dominian Astrophysical Observatory for the National Research Council of Canada's 
Herzberg Institute of Astrophysics. 
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Table 1. Submillimeter Photometry of I16293E and Total Flux Densities of I16293A, B 

Derived From Maps 





I16293E 


116293AB 






Flux Density 


Integrated Flux 


Envelope 


Filter/HPBW 


(Jy beam~^) 


(Jy beam~^) 


(Jy) 


2.0 mm/27" (2) 


0.20 ±0.03 






1.3mm/19.5" (1) 


0.60 ±0.05 


7.0 ± 1.4*^ 




1.1 mm/18.5" (3) 


0.70 ±0.07 


8.6 ±0.9 




850/xm/14.5" 




21.9 ±2.2 


16.0 


850/xm/14.0" 


1.40 ± 0.02^^ 


20.4 ± 0.03^ 




800/im/16.5" (2) 


2.06 ±0.11 


24.5 ±2.5 




750//m/13.3" 




26.6 ±10.0 




450/im/18" (2) 


11.8 ±2.2 


126 ± 15 


102.7 


450/im/8.0" 


4.35 ±0.09^ 


126 ± 0.2^ 




350/im/19" (1) 


19.4 ±2.3 


175 ± 35 





Note. — Numbers in parentheses denote the number of independent observations obtained 
for each filter/aperture combination. 

^from Mezger et al. (1992). 

''Peak flux density in cleaned SCUBA maps, see text. 

'^Integrated flux from Gaussian flt to SCUBA map. Errors are la and do not account for 
systematic calibration uncertainties. 
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Table 2. Molecular Line Observations of I16293AB and I16293E 



jTdV T Flsr 

Molecule Transition (K km s~^) (K) (km s~^) (km s~^) 



I16293AB, in units of Tmb 





3-2 


33.3 






4.0 




2-1 


18.4 


6.50 


2.50 


3.76 




2-1 


9.3 


3.24 


2.86 


3.99 




7-6 


473 






3.9 




6-5 


408 






3.9 


HCO+ 


4-3 


81 






4.1 




3-2 


70 






4.0 


H13C0+ 


4 3 


7.2 






4.7 


H13C0+ 


3-2^ 


11.8 


4.87 


2.27 


4.23 




3-2^ 


-2.3 


-2.76 


0.78 


4.24 


DCO+ 


5-4 


3.1 






4.0 




3-2 


5.18 


2.56 


2.02 


4.48 


H2D+ 


lio~lii 


0.8 


1.1 


0.6 


3.9 


N2H+ 


4-3 


4.4 


2.6 


1.7 


3.7 


HDO 




6.0 


1.0 


5.9 


4.3 




loi~Ooo 


-1.4 


-2.1 


0.60 


4.2 


H2O 


llQ-loi 


2.9 






3.9 


I16293E, in imits of Tj^ 




2-1 


5.40 


3.60 


1.41 


3.80 




2-1 


1.70 


1.26 


1.27 


4.03 


H13C0+ 


3-2 


1.64 


1.25 


1.23 


4.08 


DCO+ 


3 2 


3.74 


3.14 


1.12 


3.76 


HDO 


lio~lii 


0.1 


0.15 


0.6 


3.5 



Note. — For those spectra where we only list the integrated intensity, Vlsr is the velocity of the self- 
absorption. 

"^Emission component. 

''Absorption component. 
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Table 3. Turbulent Line Width b and Inferred Abundances for I16293AB and I16293E 







I16293A, B 


I16293A,B 


116293A,B 






I16293E 


b 


Envelope 


Envelope 


Outflow 


Species 


Abundance 


(km s^^) 


Abundance 


Abundance*^ 


Abundance 




3(-8) 


0.83 


2(-7) 


6.2(-8) 






9(-9) 


0.76 


4(-8) 


1.6(-8) 




HCO+ 


I(-IO) 


0.69 


< l(-9) 


1.4(-9) 






2(-ll) 


0.50 


2(-ll) 


2.4(-ll) 




DCO+ 


5(-ll) 


0.67 


2(-ll) 


1.3(-11) 




H2O 






3(-7)*' 
4(-9)^ 




1.3(-8) 


HDO 


2(-10) 


0.20 


3(-10) 






H2D+ 






l(-12)'i 

2(-9)'= 






N2H+ 






3(-ll) 






H3O+ 






< 5(-9) 







Note. — c{-d) denotes c x 10 

^Abundance derived by Schoier et al. (2002) in a similar analysis. 
^Where > 14 K. 
^Where Tkm < 14 K. 
'iWhere Tkin > 20 K. 
^Where Tkm < 20 K. 
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Table 4. Physical Parameters of the NE-SW Outflow 





Uncorrected^ 


Corrected'' 


Parameter 


SW Blue 


NE Red 


SW Blue NE Red 


Dynamical time scale (yr) 


6.4 X 10^ 


7.6 X 10^ 


3.0 X 10^ 3.5 X 10^ 


Mass (Mq) 


0.10 


0.05 


0.33 0.18 


Momentum (Mq km s~^) 


0.60 


0.22 


2.27 1.25 


Energy (Mq km^ s~^) 


3.44 


0.78 


23.4 8.64 


Force (10~^ Mq km s"^ yr"^) 


1.91 


0.74 


6.43 4.90 



^Uncorrected for opacity and inclination. 
'^Corrected for opacity and inclination. 
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Fig. 1.— Continuum maps of IRAS 16293-2422A, B and I16293E, all shown in greyscale 
enhanced with logarithmic contours. Top: 800 nm UKT14 map. The position of the core 
I16293E and the two protostars I16293A and B are indicated. The lowest contour is at 0.1 
Jy pixel"^ (3" pixels) and peak flux density is 1.8 Jy pixel"^. Middle: 850 /im SCUBA map. 
The lowest countour is at 0.1 Jy beam^^, peak flux is 16.0 Jy beam~^. Bottom: 450 /im 
SCUBA map. The lowest contour is at 1 Jy beam^^, peak at 76.1 Jy beam"^. The excess 
noise seen in the outskirts of the map is not real but due to edge effects. The effect from 
blanked negative bolometers north of I16293A, B is evident. The HPBW is indicated for the 
850 and 450 //m maps. 
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Fig. 2. — Observed spectra toward I16293A. 
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Fig. 3. — Observed spectra toward I16293E. 
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Fig. 4. — Azimuthally averaged flux densities of I16293A and I16293E at 850 //m {solid 
curves) and 450 //m {dashed curves). 
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Fig. 5. — Blow-up of the 450 fim. SCUBA image of I16293E. The deuterium peak position 
D (white square; Lis et ah 2002) and the peak of the SiO emission El (black dot; Hirano et 
al. 2001) are indicated. 
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Fig. 6. — Isothermal fit to the continuum emission from I16293A. For the IRAS 100 jim. data 
we show the total fiux, as well as the fiux estimated to come from the dust disk alone. The 
2.75 mm (109 GHz) point from Mundy et al. (1990) is also shown. This data point also 
includes free-free emission from both A and B and was not used in the fit. 
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Fig. 7. — Observed spectra toward I16293E overlaid with best-fit fine profiles from an isother- 
mal core model with a power-law density distribution. 



-52- 



10 



m 




10 ip 
Vlsr (km s ) 



20 




2 - 



CQ 



DCO 5-4 



H^^CO^ 4-3 



fVfff^ 




I I I I I I I I I " 




10 1^ 

Vlsr (km s ) 



20 



Fig. 8. — Selection of observed spectra toward I16293A overlaid with line profiles predicted 
by the inside-out collapse model for a sound speed a=0.7 km s^^ and age t=2.5 x 10'^ yr. A 
better fit to the wings of the C^^O and HDO spectra is obtained for a sound speed a — 0.9 
km s~^ and a — 1 km s~^, respectively {dotted curves). 



-53- 




Fig. 9. — Comparison of our power-law density and temperature distributions {solid curves) 
with those of Schoier et al. (2002) (dotted) and CeccareUi et al. (2000) {dashed curves). 
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Fig. 10. — Observed H2O spectrum {thin curves) from SWAS overlayed with our collapse en- 
velope model spectrum {top panel ); the outflow model spectrum {middle); and the coadded 
infall and outflow spectra {bottom). 
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Fig. 11. — Observed CO J =6-5 and 7-6 spectra toward II 6293 A together with overlaid 
infall {dotted) and outflow {thick) model profiles for a [C0]/[H2] abundance of 10~^. 
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Fig. 12. — Contour plots of the two outflows from IRAS 16293-2422 A and B for indicated 
velocity intervals superposed on our 800 fim continuum image. 
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Fig. 13. — Velocity position contour plots of the two outflows driven by IRAS 16293-2422 A and 
B. The top figure is created by rotating the map with the position angle of the outfiow —35° and 
integrating over a 20" wide strip approximately centered on A, i.e. this velocity position diagram 
goes along the symmetry axis of the NE-SW outflow, but also includes some contribution from 
the E-W outflow driven by B. The bottom plot is derived by rotating the CO map by 96° and 
integrating over a 20" wide strip approximately centered on B. 




Fig. 14. — H2O spectrum from SWAS {thick line) superposed on the CO J— 7-6 emission 
measured with the JCMT {dotted line). The CO emision has smoothed to the same spectral 
resolution as SWAS and scaled down to match the peak of the H2O emission. 
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Fig. 15. — Spectrum of the observed HDO ground-state transition overlaid with outflow 
model spectra for [HD0]/[H2] abundances of 3 x lO"^*^ (dotted) which matches the level of 
the maximum emission, and 3 x 10~^ (thick) which fits the outer wings. 
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Fig. 16. — Model spectra obtained from the chemical calcuations for an infalling envelope, 
overlaid on the observed H2D+ lio-lii (top) and HDO Iqi-Ooo (bottom) observations. 
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Red lobe: 
tclyn=3500 yr 
M=0.18 Msun 
[H20]=1 .3{-8) 



Prestellar core: 


[C180]=3(-8) 


p=-1.0 


[C170]=9(-9) 


Td=16 K 


[HCO+]=1(-10) 


M=4.4 Msun 


[H13CO+]=2(-11) 


Rout=8000 AU 


[DC0+]=5(-11) 




[HDO]=2(-10) 



IRAS 16293-2422 




Blue lobe: 
tdyn=3000 yr 
M=0.33 Msun 
[H20]=1.3(-8) 



Flattened core: 
R=800 AU 
Td=40 K 
M=1.8 Msun 



Infalling envelope: [C180]=2(-7) 
a=0.7km/s [C170]=4(-8) 
t=6000-25,000 yr [HC0+]<1(-9) 
M=6.1 Msun [H13CO+]=2(-11) 
R=7300AU [DC0+]=2(-11) 
Td=115...12K [N2H+]=3(-11) 



[H20, T>14K]=3(-7) 

[H20, T<14 K]=4(-9) 

[HDO]=3(-10) 

[H2D+, T>20 K]=1{-12) 

[H2D+, T<20 K]=2(-9) 

[H30+]<5(-9) 



Fig. 17. — Schematic overview of the IRAS 16293-2422 region showing its distinct physical 
regions with their derived physical parameters, and molecular abundances. 



